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The dependence of the nucleon mass on the mass of the pion is studied in the framework of the 
chiral quark-soliton model. A remarkable agreement is observed with lattice data from recent full 
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soliton model as a guideline for the chiral extrapolation of lattice data are discussed. 
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I. INTRODUCTION 

The formulation of QCD on a discrete, finite Euclidean 
lattice [1] is at present the only strict and model indepen- 
dent approach allowing to solve QCD in the low energy 
regime and to study, e.g., the hadronic spectrum from 
first principles [2]. Numerical lattice QCD simulations 
face technical problems, such as discretization errors or 
finite size effects, which are attacked and minimized with 
increasing success by employing improved versions of dis- 
cretized actions, or by working on larger lattices available 
thanks to the steadily growing computer power. Still, 
present lattices are too small to accommodate the pion 
as light as it appears in nature [3—11]. 

The tool needed to extrapolate lattice data from the 
region of nowadays typically m, > 400 MeV down to the 
physical value of the pion mass is, in principle, provided 
by the chiral perturbation theory (%PT). The xPT is 
an effective but rigorous approach to the description of 
low energy phenomena of strong interactions [12-14]. It 
is based on the concept of spontaneous chiral symmetry 
breaking with the pion as the Goldstone boson which ac- 
quires a small mass only due to explicit chiral symmetry 
breaking by the small current masses of the light quarks. 
%PT allows to address such questions like, e.g., how much 
do baryon masses change if one "switches" on the masses 
of light quarks and varies their values. 

In order to extrapolate reliably lattice data by means 
of %PT it is important to ensure the convergence of the 
chiral expansion up to large values of m % . A first and 
promising matching of x?T and lattice results was re- 
ported, and it was established that the chiral expansion is 
well under control up to m 2 < 0.4 GeV 2 [15-18], see also 
[7]. More conservative estimates, however, indicate that 
the chiral expansion is reliable only up to m 2 < 0.1 GeV 2 
[19]. The progress in computing power promisses future 
lattice data at still lower pion masses and eventually at 
the physical point, which will improve the situation and 
make disappear this problem. In the meantime, however, 
it would be desirable to have a description of an inter- 
mediate region of pion masses, that would provide a safe 
overlap between the regime of the validity of xPT and 
the lattice data. 



In this situation it is interesting to consider studies in 
chiral models - in particular, if they allow to go beyond 
the range of m w where xPT is applicable. However, the 
inevitable prize to pay for the extended range of appli- 
cability compared to xPT is model-dependence, which 
introduces hardly controlable systematic uncertainties. 
Keeping this point critically in mind, such studies may 
nevertheless provide helpful insights. 

In Refs. [20-24] the concept was introduced and de- 
veloped to regularize chiral loops by means of suitable 
vertex form factors, referred to as "finite range regula- 
tors" (FRR) and intended to simulate physical effects of 
the pion cloud which has a finite range due to m, =/= 0. 
As argued in [24] , the FRR method corresponds in some 
sense to a (though model-dependent) chiral resumma- 
tion reliable up to m 2 < lGeV 2 . While being physi- 
cally intuitive and appealing, the approach was critisized 
to be unsatisfactory from a field-theoretic point of view, 
since it gives preference of one ("finite range") regular- 
ization scheme to another (e.g. "dimensional") regular- 
ization scheme [15]. 

In this work we address the question how the nucleon 
mass depends (implicitly) on the pion mass in another 
effective approach, namely in the chiral quark soliton 
model (xQSM) [25, 26]. This model was derived un- 
der certain assumptions from the instanton model of the 
QCD vacuum [27, 28], which provides a dynamical pic- 
ture of the chiral symmetry breaking mechanism [29] . In 
this model the nucleon appears as a soliton of the chiral 
pion mean field in the limit of a large number of colours 
N c . The model provides a theoretically consistent de- 
scription of numerous baryonic quantities ranging from 
static properties [30, 31] over "usual" [32] till "general- 
ized" parton distribution functions [33] , which - as far as 
these quantities are known - agree with phenomenology 
to within (10 — 30)% at the physical point. 

Here — focussing on the nucleon mass Mjy — we 
present the first study in the xQSM at non-physical pion 
masses covering the wide range < < 1500 MeV. 
We make several remarkable observations. First, stable 
soliton solutions do exist in this range of pion masses. 
Second, we demonstrate that the model correctly de- 
scribes also the heavy quark limit. In the opposite limit 
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to w — > 0, which does not commute with large-iV c limit 
[12, 36], the xQSM is known to exhibit a chiral behaviour 
and to incorporate leading non-analytic terms, which are 
at variance with the real world QCD with a finite number 
of colours N c = 3, but in agreement with its formulation 
in the limit N c — + oo [34, 35, 37]. This is consistent 
as the model is defined in this limit. Third, we show 
that the xQSM provides a satisfactory description of the 
variation of the lattice data on Mm with in the con- 
sidered range of pion masses. 

Partly we provide explanations for these observations. 
Partly, however, they shall remain puzzles to be resolved 
upon further studies in the model. 

This note is organized as follows. In Sect. II the xQSM 
is introduced and model results for Ma^to^) presented, 
which we compare to lattice QCD in Sect. III. In Sect. IV 
we discuss the limitations for using the model quantita- 
tively to extrapolate lattice data, and compare in Sect. V 
to xPT and the FRR approach. Sect. VI contains the 
conclusions. Technical details and a digression on the 
pion-nucleon sigma-term can be found in the Appendices. 



II. PION AND NUCLEON IN THE EFFECTIVE 
THEORY 

Let us consider the effective theory which was derived 
from the instanton model of the QCD vacuum [27, 28] 
and is given by the partition function [39, 40] 

Z cS = [vipVtpVU expli [d 4 x${i@-MU j5 -m)4> ) , 



(1) 

where U = exp(iT a 7r a ) denotes the SU(2) chiral pion field 
with C/ 75 = exp(i75T a 7r°), and M is the dynamical ("con- 
stituent" ) quark mass due to spontaneous breakdown of 
chiral symmetry, and to — m u = ra^ is the current quark 
mass. We neglect throughout isospin breaking effects. 

The effective theory (1) is valid at low energies below 
a scale set by the inverse of the average instanton size 
Pav ~ 600 MeV. The dynamical mass is momentum de- 
pendent, i.e. M — M(p), and goes to zero for p >• p~*. 
In practical calculations it is convenient to replace M(p) 
by a constant mass M, and to regularize the effective 
theory within an appropriate regularization scheme with 
a cutoff A cu t = 0(p~*). In the present work we use 
M = 350 MeV from the instanton vacuum model [29] . 

In the effective theory (1) chiral symmetry is sponta- 
neously broken and a non-zero quark- vacuum condensate 
(V'V) = ( vac l(V'«V'« + VWd^vac) appears which is given 
in leading order of the large- N c limit by the quadratically 
UV-divergent Euclidean loop integral 



d 4 p E 8N C M' 



(2tt) 4 pl + M 



/2 



= -W c M'h{m), 



(2) 

where I\ is its proper-time regularized version, see 
App. A, and M' = M + m. Note that in QCD strictly 



speaking (tptp) is well-defined only in the chiral limit. The 
pion is not a dynamical degree of freedom in the theory 
(1). Instead the dynamics of the pion field appears only 
after integrating out the quark fields which yields the 
effective action 



S eS [U] = If Jd 4 x tr &>\J d^W + ... 



(3) 



where the dots denote the four-derivative Gasser- 
Leutwyler terms with correct coefficients, the Wess- 
Zumino term, terms cx m (and an infinite series of 
higher-derivative terms) [26]. The pion decay constant 
f n = 93 MeV in Eq. (3) is given in the effective theory 
by the logarithmically UV-divergent loop integral (whose 
regularized version we denote by I2, see App. A) 



r 



d 4 p E 4A C M' 2 



(2tt) 4 (p 2 e + M> 



2^2 



= 8A C M' 2 J 2 (m) . (4) 



The mass of the pion can be determined from the position 
of the pole of the pion propagator in the effective theory 
(1). Its relation to the current quark mass is given by 
the equation (for the Ii see App. A) 



2 m I\{m) 
m = — — - — - 

w M J 2 (m) 



(5) 



which, for small current quark masses to, corresponds to 
the Gcll-Mann-Oakcs-Renncr relation 



m 2 JZ = -m(W) + 0(m 2 ) . 



(6) 



The xQSM is an application of the effective theory (1) 
to the description of baryons [25, 26]. The large- N c limit 
allows to solve the path integral over pion field configura- 
tions in Eq. (1) in the saddle-point approximation. In the 
leading order of the large- N c limit the pion field is static, 
and one can determine the spectrum of the one-particle 
Hamiltonian of the effective theory (1) 

H\n) = E n \n) , H = -ij°j k d k + 7 °MC/ 75 + 7 °to . (7) 

The spectrum consists of an upper and a lower Dirac 
continuum, distorted by the pion field as compared to 
continua of the free Dirac-Hamiltonian 

H \n ) = E no \n ) , H = -i 7 ° 7 fe a fe + 7 M + 7 to , (8) 

and of a discrete bound state level of energy E\ ev , if the 
pion field is strong enough. By occupying the discrete 
level and the states of the lower continuum each by A^ c 
quarks in an anti-symmetric colour state, one obtains a 
state with unity baryon number. The soliton energy E so \ 
is a functional of the pion field 



E sol [U] = N c 



E\ev + ^ {E n - E na ) 

E n <0 



(9) 



Ssoi[f] is logarithmically divergent, see App. A for the 
explicit expression in the proper-time regularization. By 
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minimizing _E so i[£7] one obtains the self-consistent soli- 
tonic pion field U c . This procedure is performed for sym- 
metry reasons in the so-called hedgehog ansatz 

7r a (x) = e% P{r) , C/(x) = cos P{r) +ir a e a r sin P(r) , 

(10) 

with the radial (soliton profile) function P(r) and r = |x|, 
e r = x/r. The nucleon mass Mat is given by E so i[U c ]. 
The self-consistent profile satisfies P c (0) = — tt and de- 
cays in the chiral limit as 1/r 2 at large r. For finite m 
it exhibits a Yukawa tail oc cxp(—m n r)/r with the pion 
mass TOtt connected to m by the relation (5). 

For the following discussion we note that the soliton 
energy can be rewritten as an expansion in powers of 
d^U as follows 



E sol [U) = Y^F k [{dUf 



(11) 



fc=i 



where Fk[(dU) k ] symbolically denotes a functional in 
which d^U appears fc-timcs (appropriately contracted). 
Note that in leading order of the large- iV c limit the soli- 
ton field is time- independent, i.e. dU is just VU. For 
some observables the lowest orders in expansions analog 
to (11) were computed [30-33]. 

Let us also remark that in the case m u = m d = m s the 
above given formulae for the soliton energy in the SU(2) 
version of the model coincide with those from the SU(3) 
version [31]. 

In order to describe further properties of the nucleon, 
it is necessary to integrate over the zero modes of the 
soliton solution in the path integral, which assigns a defi- 
nite momentum, and spin and isospin quantum numbers 
to the baryon. Corrections in the 1 /./V c -expansion can be 
included by considering time dependent pion field fluctu- 
ations around the solitonic solution. The xQSM allows 
to evaluate baryon matrix elements of local and non-local 
QCD quark bilinear operators like (S'|*r*|B) (with T 
denoting some Dirac- and flavour-matrix) with no ad- 
justable parameters. This provides the basis for the wide 
range of applicability of this model [30-33] . 



Mjv ~ N c Mi = E so i. We obtain numerically 1 
M N = 1254 MeV 



(13) 



where the cutoff in the regulator function R, see App. A, 
is chosen such that for m v = 140 MeV the physical value 
of the pion decay constant is reproduced. We observe 
an overestimate of the physical value = 940 MeV 
by about 30%. This is not surprising, given the fact 
that we truncate the series in Eq. (12) after the first 
term and thus neglect corrections which are generically 
of 0(1/N C ), i.e. of the order of magnitude of the observed 
overestimate. 

In fact, the soliton approach is known to overestimate 
systematically the physical values of the baryon masses 
because of - among others - spurious contributions due 
to the soliton center of mass motion [38] . Taking into ac- 
count the corresponding corrections, which are <D(N®), 
reduces the soliton energy by the right amount of about 
300 MeV. (Note that there are also other sources of cor- 
rections at 0(N°), see Ref. [38].) We shall keep in mind 
this systematic overestimate, when we will discuss lattice 
data below. 



B. The chiral limit 

In the following we will be interested in particular in 
the pion mass dependence of the nucleon mass. From 
xPT we know that 



M N (ra,) = M N (0) + A ml + B ml + . . . 



(14) 



where the dots denote terms which vanish faster than 
ml with TOtt — > 0. The constants Mn(0) and A (which is 
related to the pion-nucleon sigma-term) serve to absorb 
infinite counter-terms in the process of renormalization in 
xPT - in the sense of renormalizability in xPT. However, 
the constant B, which accompanies the so-called leading 



non-analytic (in m, since ml 



,3/2 



) contribution, is 



finite. For this constant xPT, as well as any theory which 
correctly incorporates chiral symmetry, yields [12, 36] 



A. The mass of the nucleon in the large-iV c limit 

If one managed to solve QCD in the limit N c — > oo one 
would in principle obtain for the mass of the nucleon an 
expression of the form (let the Mi be independent of N c ) 



M N = N C M\ + N°M 2 + N^Ms + ... 



(12) 



B 



k 



3gj 

327T/2 



(15) 



with k = 1 for finite N c . However, the limits m^ — > 
and N c — > oo do not commute. If we choose first to 
take N c — > oo while keeping m T finite, and only then 
we consider the chiral limit, then k = 3. The reason for 
that is the special role played by the A-resonance. In the 



The xQSM provides a practical realization of the large- 
N c picture of the nucleon [41], and respects the general 
large- N c counting rules. In leading order of the large- N c 
limit one approximates in the xQSM the nucleon mass 
by the expression for the soliton energy in Eq. (9), i.e. 
one considers only the first term in the expansion (12) 



1 In this work we quote numerical results to within an accuracy of 
1 MeV. However, one should keep in mind that we neglect isospin 
breaking effects (and electromagnetic corrections). Therefore we 
shall round off the physical masses as = 140 MeV and Mjv = 
940 MeV, and the same is understood for our numerical results. 
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large- N c limit the nucleon and the A-resonance become 
mass-degenerated 2 



Ma — Mn — OiN- 1 ) . 



(16) 



Taking N c — > oo while m T is kept finite, one has 
Ma — Mjy <C m T and has to consider the contribution of 
the A-resonance as intermediate state in chiral loops on 
equal footing with the contribution of the nucleon. The 
contribution of the A-resonance appears (in quantities 
which do not involve polarization) to be twice the con- 
tribution of the nucleon, hence k = 3 [37]. This is the 
situation in the xQSM and, in fact, in this model one 
recovers [35] the correct (in the large-7V c limit) leading 
non-analytic term in Eq. (15). 



C. The heavy quark limit 

The xQSM is based on chiral symmetry which consid- 
eration makes sense only when explicit chiral symmetry 
breaking effects are small. Nevertheless one can consider 
the model, in principle, for any value of m. In this context 
let us first note that taking the limit of a large current 
quark mass m — ► vhq in Eqs. (2, 4, 5) yields, see Table I 



lim m„ 

m — >m,Q 



2m Q 



(17) 



which is the correct heavy quark limit for the mass of a 
meson (see Sec. II D on details of parameter fixing). 

What happens to the mass of the soliton in this limit? 
Consider the eigenvalue problem in Eq. (7). With in- 
creasing m^> niQ the "potential term" MJ7 75 is less and 
less important and the spectrum of the full Hamiltonian 
(7) becomes more and more similar to the spectrum of the 
free Hamiltonian (8). In this limit the eigenvalues of the 
full and free Hamiltonian nearly cancel in the sum over 
energies in Eq. (9). Only the contribution of the discrete 
level is not compensated and approaches E\ ev w tuq. 
Thus, we recover the correct heavy quark limit 



lim Mn = N c mn . 



(18) 



However, the above considerations are formal, since the 
crucial step consists in demonstrating that a stable soli- 
ton solution indeed exists, i.e. that a self-consistent pro- 
file can be found for which the soliton energy (9) takes 
a minimum. In practice, we find that stable soliton so- 
lutions exist at least up to m = 0(700 Me V), which is 



2 Though in the xQSM one cannot handle 0(N®) corrections to 
Mat, still one is able to describe consistently the mass differ- 
ence (Ma — Mpf) and to reproduce the large-iV c counting rule 
in Eq. (16) by considering a class of particular (so-called "ro- 
tational") l/7V c -corrections [28]. In the soliton picture A and 
nucleon are just different rotational states of the same classical 
soliton. In this note we do not consider rotational corrections 
and work to leading-order in the large- N c expansion. 



M N ( m7l )/GeV 



exact result, here 
approximation ofRef.f 14] 



m^GeV 2 



FIG. 1: Nucleon mass Mpf(m 7v ) vs. m 2 in the chiral quark-soliton 
model. Solid line: Exact result obtained here. Dashed line: The ap- 
proximation based on the instanton vacuum model from Ref. [35]. 



sufficient for a comparison to lattice QCD results. In 
this range we observe a tendency to approach the limit 
(18) from below as expected, see Table I. Thus, we find 
that both the pion and the nucleon mass are correctly 
described in the effective theory when the current mass 
of the quarks becomes large. 



D. Fixing of parameters 

When studying the nucleon mass as function of we 
must specify what is fixed and what varies in the chiral 
limit. Here we make the following choice. We keep the 
dynamical constituent quark mass M = 350 MeV and 
/tt = 93 MeV fixed. In this way we obtain the results 
shown in Table I, and plotted as solid line in Fig. 1. The 
meaning of the dashed curve in Fig. 1 is explained in 
App. B where a digression is given on the pion- nucleon 
sigma-term related to the slope of Mjv(m w ). 

It should be noted that keeping f n fixed in the chiral 
limit is a choice often considered in literature. However, 
at this point other choices could be considered as well. 
E.g. one could fix the pion decay constant to its value 
F = 88 MeV in the chiral limit, or allow f n to be m n - 
dependent, which strictly speaking is the case in xPT 
and in lattice QCD. In xPT 



U{m*) = F 1 + 



(4ttF)< 



(19) 



where I4 is a low energy constant [13, 14]. In lattice 
calculations f n increases with larger m v , exceeding its 
physical value by about 40% at ~ 1 GeV, e.g. [42]. 
A more consistent way of fixing model parameters could 
consists in choosing A cut (and/or M) such that in the 
model f n (m n ) satisfies (19) with the correct value for I4, 
and agrees in each case with lattice results at large m,. 

Remarkably, Eq. (19) holds in the model and parame- 
ters can be fixed to reproduce I4 correctly [34]. However, 
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it is a subtle issue, how to simulate in a chiral model 
the lattice situation in a realistic way. There a common 
procedure is to keep fixed all bare lattice parameters but 
the bare current quark mass (the hopping parameter) . In 
some sense, the dimensionful quantity, which is kept fixed 
in lattice calculations while m n "is varied", is the Som- 
mer scale [43] defined by the notion of the heavy quark 
potential — absent in chiral models. 

Notice that the xQSM describes numerous observables 
to within an accuracy of typically (10 — 30)% [31], but 
cannot be expected to be much more precise than that. 
Overestimating the value of f w in the chiral limit by 5% 
or underestimating its lattice values by 30% (or 40%) lies 
within the accuracy of the model. Note that changing 
(at non-physical m^) would alter in particular values of 
the cutoff A cut and consequently change in Table I. 
Since and Mn depend on A cut logarithmically and 
thus weakly, these changes may be expected to be small 
and within the accuracy of the model. 

From a practical point of view, our choice to keep 
fixed to its physical value may be considered as one ef- 
fective prescription, whose consequences one may think 
absorbed in the unavoidable model dependence of the 
results. 

It would be interesting to consider other effective pre- 
scriptions which would provide more insights into the 
model dependence of our study. This is, however, beyond 
the scope of this work and subject to further investiga- 
tions [44]. 



TABLE I: The dependence on the pion mass for fixed = 93 McV 
in the xQSM. All numbers are in units of McV. Rows 2,3,4: Current 
quark mass m, cutoff of the effective theory A cu t, and quark vac- 
uum condensate (^ij)) depend on according to Eqs. (2, 4, 5). 
A cu t is of the order of magnitude of the inverse of the average 
instanton size p^ v 600 Me V. JVote that in QCD - in con- 
trast to effective theories with a well defined regularization pre- 
scription - the notion of quark vacuum condensate for m ^ 
is ambiguous. Rows 5,6,7: Contributions of the discrete level 
N c E\ ev and the continuum N c E con t to the total soliton energy 
Esol = N c (Ei cv + -Econt); see Eq. (9), to be identified with the nu- 
cleon mass in the model. The numerical numbers confirm within 
the studied range of m the heavy quark limit discussed in Sec. UC. 
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III. COMPARISON TO LATTICE DATA 

In this Section we compare MN(m 7: ) from xQSM with 
lattice results on the nucleon mass from Refs. [4-11]. 
Wherever possible we will use only such lattice data, 
where all bare lattice parameters where kept fixed apart 
from the current quark mass (or the hopping paramater). 

M]v(m T ) from the model overestimates the lattice data 
by about 300 MeV which is expected, if we recall the 
discussion in Sec. II A. Remarkably we observe that it 
is possible to introduce an m^-independent subtraction 
constant C — O(300 MeV) depending on the lattice data, 
such that [Mxim^) — C] agrees well with the lattice data. 

In Fig. 2a we compare the xQSM result for Ma^to^) to 
full lattice data by the UKQCD and QCDSF Collabora- 
tions [6, 7] obtained from simulations on (16 — 24) 3 x 48 
lattices using the standard Wilson plaquette action for 
gauge fields and the non-perturbatively 0(a) improved 
action for fcrmions. The fat points in Fig. 2a were ex- 
tracted from simulations at 5.20 < [3 < 5.29 on lattices 
of the physical size L > 1.96 fm with m^L > 6 and lat- 
tice spacings a — 0.09 fm to 0.12 fm fixed by the Som- 
mer method with r = 0.5 fm [43]. Other lattice spacings 
quoted below were also determined by means of this pop- 
ular method which, however, is not free of criticism [45] . 
A best fit yields C = 329 MeV with a x 2 per degree of 
freedom of xlof = 0-7 and is shown in Fig. 2a. The un- 
certainty of the constant C due to the statistical error 
of the lattice data is of the order of magnitude of few 
MeV, see Table II and thus negligible, see the remarks in 
footnote 1. 

For comparison in Fig. 2a also data by the UKQCD 
and QCDSF Collaborations [6, 7] are shown from smaller 
lattices L = 1.5 fm to 1.7 fm. Finite size effects are clearly 
visible, and were subject to a detailed study in Ref. [7]. 
Since the present xQSM calculation by no means is able 
to simulate finite volume effects, we restrict our study to 
data from lattices with L > 2 fm. The study of Ref. [7] 
indicates that for lattices of this size finite volume ef- 
fects can be assumed to be small. Also, in this work 
we consider data obtained from lattices with spacings 
a < 0.13 fm and assume discretization effects to be negli- 
gible. Such effects are difficult to control systematically 
[46]. 

In Fig. 2b the xQSM result for M N (m^) is com- 
pared to lattice data by the JLQCD Collaboration from 
Ref. [5], where dynamical simulations with two flavours 
were performed on a 20 3 x 48 lattice using the plaque- 
tte gauge action and the non-perturbatively 0(a) im- 
proved Wilson quark action at (3 — 5.20 with lattice 
spacings 0.10 fm < a < 0.13 fm. The physical size of 
the lattice was 1.96 fm < L < 2.53 fm and the range 
500 MeV < m n < 1 GeV was covered corresponding to 
m^L > 5. The best fit yields C = 341 MeV with a 
Xdof = 0-3 and is shown in Fig. 2b. 

In Figs. 2c and d we compare M N {m^) from the 
model to lattice results by the CP-PACS Collabora- 
tion [4], which were obtained on 24 3 x 48 lattices from 
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(a) 



QCDSF-UKQCD Lattice Collaborations 
%QSM: Mf/mJ - 0.328 GeV 



- extrapolated M N = 925 MeV 



2 m^/GeV 2 



M N ( mjt )/GeV 



(c) 



CP-PACS Collaboration: p = 2.1 
%QSM: Mf/mJ - 0.331 GeV 



- extrapolated M N = 922 MeV 



2 m^/GeV 2 



M N (m^)/GeV 
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2 - 



(b) 



JLQCD Collaboration 
%QSM: Mf/mJ - 0.341 GeV 



- extrapolated M N = 913 MeV 



2 m^GeV 2 



M N (m^)/GeV 



(d) 



CP-PACS Collaboration: (, - 2.2 
%QSM: Mf/mJ - 0.338 GeV 



- extrapolated M N = 915 Me V 



2 m^GeV 2 



M^m^/GeV 
3 



(e) 



RBC Collaboration 
%QSM: M^/mJ - 0.306 GeV 



- extrapolated M N = 947 MeV 



2 m^GeV 2 



M^m^/GeV 
3 



(f) 



MILC Collaboration, full (2001 ) 
MILC Collaboration, full (2004) 
XQSM: Mf/mJ - 0.256 GeV 



~ extrapolated M N = 998 Me V 



2 m^GeV 2 



FIG. 2: Mf](m- K ) vs. from the chiral quark-soliton model vs. lattice data on Mjv(nijr) from full simulations by (2a) the UKQCD- 
QCDSF Collaborations [6, 7], (2b) the JLQCD Collaboration [5], (2c,d) the CP-PACS Collaboration [4], (2e) the RBC Collaboration [8], 
(2f) the MILC Collaboration [9, 10]. All data were obtained from large lattices L > 2fm, with the exception of the data marked by crosses 
in (2a), see text. Unless error bars are shown, here and in the following figures the statistical error of the lattice data is comparable to or 
smaller than the size of the points in the plot. 




FIG. 3: Mjv(m^) vs. m\ from the chiral quark-soliton model vs. lattice data on Mj^(m^) obtained in the quenched approximation by 
3a the CSSM Lattice Collaboration [11] and 3b the MILC Collaboration [9]. 
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a renormalization-group improved gauge action and a 
mean field improved clover quark action at f3 = 2.1 (and 
P = 2.2) with dynamical u- and d-quarks and quenched 
s-quarks. For hadrons with no valence s-quarks this prac- 
tically means a full two-flavour simulation. The physical 
lattice spacings and sizes were a = 0.09 fm to 0.13 fm 
and L — 2.22 fm to 3.12 fm with pion masses in the 
range 500 MeV < m v < 1 GeV, such that m^L > 7. 
We observe that Mjv(m w ) from the xQSM with a sub- 
traction constant C = 331 MeV with x\ a f = 2-2 (and 
C = 338 MeV with xlof = °- 4 ) describes well the lattice 
data of Rcf. [4], see Fig. 2c (and Fig. 2d). 

Next we confront the model results for M N (m^) to 
lattice data from dynamical two-flavour simulations with 
domain wall fermions by the RBC Collaboration [8], 
which have the virtue of preserving chiral invariance. 
In Ref. [8] a renormalization group improved ("doubly 
blocked Wilson") gauge action with (3 = 0.80 was used 
on a 16 3 x 32 lattice with the physical spacing of 0.12 fm 
and a lattice size about 2 fm. The range of pion masses 
was win = 470 MeV to 650 MeV. A best fit yields for the 
constant C = 306 MeV with a x\ f = 0-02, and provides 
a very good description of the data, see Fig. 2e. 

In Fig. 2f we confront M^TOtt) from the xQSM with 
2001 lattice data by the MILC Collaboration [9] ob- 
tained from simulations with three dynamical quarks 
using a one-loop Symanzik improved gauge action and 
an improved Kogut-Susskind quark action. The phys- 
ical lattice size was tuned to L = 2.6 fm with a lat- 
tice spacing of 0.13 fm and the range of pion masses 
340 MeV < nix < 2.2 GeV was covered. The best fit 
yields C = 256 MeV for the subtraction constant (with 
XL = 1-2). 

One could worry whether the SU(2) model results 
can be compared to three-flavour lattice simulations, 
though in [9] for the nucleon mass no significant differ- 
ences were noticed between two- and three-flavour runs. 
However, as noted in Sec. II, the nucleon mass is the 
same in the SU(2) and SU(3) versions of the xQSM if 
m u = ma = m s , which is the case for the lattice data [9] 
for TOtt > 700 MeV. Restricting the fitting procedure to 
this range of m w (the last three points in Fig. 2f) we ob- 
tain C = 260 MeV (with a x\ oi = 1-9)- This fit is shown 
in Fig. 2f. Note that it is practically indistinguishable, 
c.f. footnote 1, from the fit where the whole m^-range (i.e. 
also data with m u — < m s ) was used, and equally 
well describes the region of lower m v . In any case we 
observe a good agreement with the lattice data [9] up to 
m v = 1.5 GeV. 

We include in Fig. 2f also the recent small-m^ (2004-) 
MILC data [10] from the "coarse" lattices with L = 20 
(or L = 24 for the lowest m^-value) which have a lattice 
spacing a w 0.12 fm comparable to [9]. These data are 
not compatible with the xQSM result. In this context it 
is important to note that the simulation for the highest 
m^-value of the 2004 MILC data [10] is an extended run 
of the simulation for the lowest of the 2001 data [9]. 
One would therefore expect that they coincide in the plot 



in Fig. 2f, which is not the case. In fact, the MILC data 
[9, 10] for this particular simulation are well consistent 
with each other in lattice units: am^ and aMjy from 
these runs agree within statistical error bars. 

The discrepancy of the results from [9, 10] in Fig. 2f 
is due to the different values for r\ — a parameter de- 
fined similarly to the Sommer scale ro and used to fix the 
physical units in [9, 10]. Different methods were used to 
determine the physical unit of ri — resulting in the value 
n = 0.35 fm in [9], and n = 0.324(4) fm in [10]. Had we 
used the ri-value from [9] to give physical units to the di- 
mcnsionless lattice numbers for am n and aMjy from [10], 
then the two data sets would perfectly agree in Fig. 2f, 
and the 2004 data would be compatible with the xQSM 
results. The precise determination of the physical units 
of lattice data is a difficult issue [43, 45], see also [9, 10]. 

It is instructive to compare the model results also 
to lattice data obtained from simulations performed in 
the quenched approximation, e.g., by the CSSM Lat- 
tice Collaboration [11], where a mean-field improved 
gauge action and a fat-link clover fermion ( "FLIC" ) ac- 
tion was used on a 16 3 x 32 lattice with a lattice spac- 
ing of a = 0.125 fm. The calculation covers the range 
540 MeV < < 920 MeV. A best fit to the quenched 
data gives C = 345 MeV (with a Xdof = 1-3) which yields 
a good agreement with the lattice data, see Fig. 3a. 

For the quenched data by the MILC Collaboration [9] , 
however, we observe that a fit would work much worse. 
In this case we refrain from fitting and show instead in 
Fig. 3b the fit to the full MILC data from Fig. 2f, which 
nicely illustrates how results from full and quenched cal- 
culations differ. Interestingly, at large the full and 
quenched data of Ref. [9] agree well with eachother. In 
fact, it is not suprising that differences between full and 
quenched simulations become less pronounced with in- 
creasing m^., i.e. with increasing fermion masses. 

Thus, we observe that the xQSM-result for Ma^to^) 
supplemented by an m^-independent subtraction con- 
stant C (whose precise value follows from a best fit to 
the respective lattice data) is able to describe the lat- 
tice data [4-9, 11] over a wide range of pion masses 
350 MeV < TOtt < 1500 MeV. The results are shown in 
Figs. 2 and 3 and are summarized in Table II. 



IV. xQSM AS TOOL FOR EXTRAPOLATION? 

From the [M N (m n ) - C] at m n = 140 MeV with the 
constant C fitted to the respective data, we can read 
off in principle the physical value of the nucleon mass - 
"extrapolated" from the respective lattice data by means 
of the xQSM as a guideline. These extrapolated values 
arc indicated by arrows in Figs. 2a-f and 3a, and arc 
summarized in Table II. 

It is worthwhile stressing that the extrapolated values 
in Table II agree within an accuracy of ±5% with the 
physical nucleon mass. (Though at the same time, lat- 
tice data for / w are underestimated by up to 40% at large 
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TABLE II: Comparison of the M^(m n ) obtained from the xQSM to the lattice data [4-9, 11]. For convenience we quote the lattice sizes 
and spacings and the range of covered in the lattice simulations in physical units which were fixed by the Sommer method [43]. The 
soliton approach generally overestimates [38] the nucleon mass at the physical point by about 300 MeV, see Eq. (13). We find a similar 
overestimate at the respective lattice values of m w . Correcting for this overestimate by introducing a m^-independent subtraction constant 
C to be fitted to the respective lattice data set, we observe a good agreement [Mjv(m.7r) — C] with the lattice data, see Figs. 2 and 3. The 
5th row shows the fit results for the constant C and its 1-rr uncertainty due to the statistical error of the lattice data, and the 6th row 
shows the x 2 per degree of freedom (x\ Q f) of the respective fit. Also the "extrapolated" (within the xQSM) value of the nucleon mass at 
the physical point (Mn ) is included. It has the same uncertainty as the fit-constant C, which is due to the statistical error of lattice data 
and practically negligible, see the remark in footnote 1. It has also an uncstimatcd systematic error due to model-dependence, see Sec. IV. 



Collaboration 


lattice size/fm 


spacing/fm 


m^-rangc/MeV 


C/MeV 


2 

Xdof 


M w /MeV 


UKQCDSF [6, 7] 


2.0-2.2 


0.09-0.12 


550-940 


329 ± 7 


0.7 


925 


JLQCD [5] 


2.0-2.5 


0.10-0.13 


540-960 


341 ± 7 


0.3 


913 


CP-PACS [4] (/3 = 2.1) 


2.7-3.1 


0.11-0.13 


520-960 


331 ± 4 


2.2 


922 


CP-PACS [4] (13 = 2.2) 


2.2-2.4 


0.09-0.10 


590-970 


338 ± 6 


0.4 


915 


RBC [8] 


1.9-2.0 


0.12 


470-650 


306 ± 11 


0.02 


947 


MILC [9] (all data) 


2.6 


0.13 


350-1500 


256 ± 2 


1.2 


998 


MILC [9] (only m u , d = m a ) 


2.6 


0.13 


720-1500 


260 ± 3 


1.9 


993 


CSSM [11] (quenched) 


2.0 


0.125 


540-920 


345 ± 5 


1.3 


909 



TOtt <~ 1 GeV, see Sec. II D.) The uncertainty of the ex- 
trapolated values for the nucleon mass is the same as for 
the fit constant C as quoted in Table II, i.e. it is of the 
order of few MeV and thus negligible, see footnote 1. 

However, what does not need to be negligible, is the 
systematic error of such an extrapolation. First of all, 
this extrapolation is done within the xQSM, and the re- 
sults are model-dependent. E.g. handling the model pa- 
rameters in the chiral limit differently would change our 
results — though one may expect a qualitatively simi- 
lar picture, see the discussion in Sec. II D. Apart from 
this source of model-dependence, which will be subject 
to future numerical studies [44], there are principal dif- 
ficulties to estimate reliably the systematic error within 
the model, as we shall see in the following. 

Let us first address the question of the range of relia- 
bility of the xQSM-description of Mjv(m,r). For that it 
is instructive to compare to x^T and the effective FRR 
approach. xPT to 0(p 4 ) was argued to provide a reli- 
able expansion for Mj^(m n ) up to m 2 < 0.4 GeV [15-18] 
(p is to be identified with the generic small expansion 
parameter in xPT, e.g., in our context pion mass). A 
more conservative bound m 2 < 0.1 GeV 2 was given in 
[19]. The effective FRR approach was argued to corre- 
spond to a partial resummation of the chiral expansion, 
and to be valid up to m 2 < 1 GeV 2 [24] (which comes, 
of course, at the prize of introducing model dependence, 
see the discussion in Sec. I.) 

What might be the range of reliability of the xQSM? 
Recall that the xQSM expression for the nucleon mass 
(as well as for any other quantity) may be considered 
as a resummed infinite series in derivatives of the chiral 
field U = exp(ir a 7r a ), c.f. Eq. (11), whereby it is under- 
stood that each chiral order is evaluated in the large- N c 
limit. (Keep in mind that these limits do not commute, 



see Sect. II B.) Thus, one may argue that the xQSM 
corresponds to a chiral expansion, which is completely 
resummed - to leading order of the large- N c expansion. 
If one were happy with this approximation, then xQSM 
results could be considered reliable for all m T including 
the heavy quark limit, see the discussion in Sec. II C. It 
must be stressed, that - as in the case of the FRR ap- 
proach - this chiral resummation is performed within a 
particular model of the nucleon. Thus, our results and 
conclusions are inevitably model dependent. 

The wide range of reliability of the xQSM we observe 
in practice is illustrated in Fig. 4, where we show the 
xQSM-fit to the MILC lattice data [9] and the xPT-fit 
from [17]. (For further details see Sec. V.) 

The perhaps most serious restriction might be that the 
xQSM assumes the number of colours N c to be large. 
What might be the effects due to l/V c -corrections? 

M N (m n )/GeV 

I , , , , | , , , , | , , , , | I 

^ _ MILC 

XQSM: MrfmJ- 0.256 GeV 

%PT to 0(p") from Ref.[16] 



XQSM-extrapolated 
M N - 998 MeV 




mJGeV 



FIG. 4: Lattice data by MILC [9] on M^r(m^) vs. and the fits 
to these data in xPT horn Ref. [17] (where the physical value of 
Mjy was used as input) and in the xQSM, c.f. Fig. 2f. The figure 
illustrates the wide range of applicability of the xQSM. 
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Let us consider the mass difference between the A- 
resonance and the nucleon as a measure for such correc- 
tions. Note that this mass difference vanishes not only in 
the large- N c limit, see Eq. (16), but also in the heavy 
quark limit [36]. This is supported by lattice results 
[4, 9, 11], where always A = Ma - Mn < m n holds. 
(This means that on present day lattices the A-resonancc 
is safe from strong decays and thus a stable particle.) In 
the region m 2 > 0.32 GeV 2 one finds A < \m v [4, 9, 11]. 
Thus, practically for most of the present day lattice data 
the condition 

A = Ma — Mn -C for to 2 > 0.32 GeV 2 (20) 

is satisfied, such that the A-nucleon mass difference can 
be neglected to a good approximation. This may be a rea- 
son for the good description of lattice data in the xQSM 
up to ml < 2.3 GeV 2 in Sec. III. 

However, in the physical region A = 290 MeV is larger 
than m^ = 140 MeV, and the leading order large- N c 
treatment of the nucleon mass inevitably introduces a 
serious systematic error. In fact, one could express the 
nucleon mass as a function of m^ and y = A/m, as 

M N (m w ,A) =F(m v ,y) , y= — . (21) 

Then, using the xQSM as a guideline for the extrapola- 
tion of lattice data corresponds to approximating 

M N = F(m w ,2.1) w F(m n ,0) at m^ = 140 MeV. 

(22) 

It is difficult to quantify the systematic error associated 
with this approximation. A very rough estimate of this 
error within the model is given in Appendix C. 

Finally, let us discuss the role of the subtraction con- 
stant C. As mentioned in Sec. II the appearence of 
such a constant is theoretically well motivated and un- 
derstood in the soliton approach [38]. The compari- 
son to lattice data indicates that this constant is about 
C = O(300 MeV) and m^-independent in the covered 
m^-range within the statistical accuracy of the lattice 
data, see Table II. Although this happens to be the mag- 
nitude for this constant needed for the model result to 
coincide with the physical value of the nucleon mass, see 
Eq. (13), it would be premature to assume the constant C 
to be m^-independent for all m v . However, this is what 
we implicitly we did when quoting the extrapolated val- 
ues for M/v in Table II. The question, whether or not 
the constant C is m n dependent, cannot be answered 
rigorously within the model. 

As an intermediate summary, we conclude that besides 
the general drawback of being model-dependent, the use 
of the xQSM as a guideline for an extrapolation of lattice 
data is limited by two major sources of systematic error, 
namely 1/N C corrections and a possible m^-dependence 
of the constant C . Both are not under control within the 
model and prevent a reliable estimate of the systematic 
error of the extrapolation. At this point it is instructive 
to compare to xPT and the FRR model, which may give 
us a rough idea about the size of the systematic effects. 



V. COMPARISON TO X PT AND FRR 

In order to test the results from the xQSM below the 
m n range available from lattice QCD we have to compare 
to xPT which allows to connect model-independently lat- 
tice data through the physical point to the chiral limit. 
In the following we shall focus on the analyses [7, 16, 17] 
in xPT up to 0(p ) (and refer to them simply as xPT). 
To this order the chiral expansion of the nucleon mass is 
characterized by 4 low-energy constants — or more pre- 
cisely, by 4 linearly independent combinations of them. 
These constants are in principle known from studies of 
nucleon-nucleon or pion-nucleon low-energy scattering 
data. However, they can alternatively be determined 
from a fit to lattice data. 

Fig. 5a shows the Mjv(m,r) as obtained in x?T from 
fits to the lattice data [4-7] satisfying a < 0.15 fm, 
m^L > 5 and respectively the constraint m n < 800 MeV 
in Ref. [7], and m T < 600 MeV in Rcf. [16]. Also shown 
in Fig. 5a is the fit of Ref. [17] to the MILC data [9], 
where a simultaneous fit to lattice data on Mi(m,) was 
performed. In [16, 17] the physical values of the nucleon 
and/or A- mass were included as constraints to the fit. 
Note that low-energy constants resulting from the dif- 
ferent fits [7, 16, 17] arc compatible with the respective 
phenomcnological values. From this point of view the 
different lattice data in Fig. 5a are consistent with each 
other. For comparison in Fig. 5a our xQSM-fits to the 
same data sets are shown. 

Notice that the two highest-m w data points (marked 
by empty circles in Fig. 5a) by the UKQCD Collabora- 
tion [6, 7] were obtained from somehow smaller lattices 
of the size L = 1.56 fm and 1.68 fm and clearly do not fol- 
low the tendency of the M^m^) from the xQSM. These 
points are therefore omitted from the fit shown in Fig. 5a. 
Including these points significantly worsens the x 2 f of 
the fit, see Table III where the results are summarized. 
Since the xQSM-description of the lattice data effectively 
works also at significantly larger m^, see Sees. Ill and 
IV and c.f. Fig. 4, we conclude that, quoting Ref. [16], 
"the surprizingly good (and not yet understood) agree- 
ment with lattice data [above m w > 600 MeV] even up to 
m n ~ 750 MeV" is an accidental consequence of compar- 
ing xPT to 0(p A ) at the edge (if not above) its range of 
applicability to lattice data where finite size effects start 
to play a role. 

In fact, up to m v < (500 - 550) MeV the xPT to 
0(p A ) [7, 16, 17] and the xQSM describe M N (m n ) in 
good qualitative agreement, see Fig. 5a. Beyond this 
point, however, the nucleon mass as function of m 2 from 
xPT in Refs. [7, 16] changes the curvature, indicating 
that the range of reliability of x?T to 0(p 4 ) could be 
m 7r ^ 500 MeV (which, in fact, is not far from the gener- 
ally assumed bound m w < 600 MeV) . This contrasts the 
xQSM-result exhibiting in agreement with lattice data 
negative curvature up to the highest considered to 2 . 

With the above considerations in mind, we conclude 
that the xQSM describes lattice data [4-7] constrained 
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FIG. 5: (a) Selected lattice data on M N (m n ) by QCDSF- UKQCD, CP-PACS, JLQCD [4-7] and the corresponding fits in X PT to 0(p 4 ) 
from Rcfs. [7, 16], in comparison lattice data by MILC [9] and the fit in \PT to 0(p 4 ) from Rcf [17]. (b) Lattice data by CP-PACS [4] 
and the fit in the finite range regulator approach (with a dipole regulator) [24]. In both figures the respective xQSM-hts are shown. 



by a < 0.15 fm, m^L > 5 and L > 2 f m very well, see 
Fig. 5a and c.f. Fig. 4. In the range below < 500 McV 
we observe a good qualitative and quantitative agreement 
of the xQSM with xPT to 0(p A ) [7, 16, 17]. In this 
range the curves for Mm(tti t ) from the two approaches 
agree with each other to within an accuracy of 50MeV 
and better, see Fig. 4 and 5a. This number may give 
us a flavour of the magnitude of the systematic error of 
the xQSM-extrapolation of lattice data, though x?T to 
0(p ) may also have such an intrinsic uncertainty, as 
indicated in Table III, due to unestimated contributions 
from 0(p 5 ). 

Next we consider the effective approach based on the 



TABLE III: The value of the physical nuclcon mass Mjv as extrap- 
olated by means of\PT to 0(p 4 ) lattice data [4-7] subject to the 
constraints a < 0.15 fm, m^L > 5 and respectively m n < 800 McV 
in Ref. [7], and < 600 MeV in Ref. [16]. For comparison the 
XQSM-hts to the same data set, and to the same data set subject 
to the additional condition L > 2 fm are shown. The error of Mjv 
due to the statistical uncertainty of the lattice data and the Xdof °^ 
the fit are shown. The systematic error of Mjv due to the extrapo- 
lation method, as well as lattice discretization or finite size effects, 
is merely indicated. (See the remark in footnote 1 concerning the 
small error of the xQSM-fit due to the statistical uncertainty of 
lattice data.) 



Method Mjv in MeV Xdof 



X PT to 0(p 4 ), "Fit I" in Rcf. [7] 948 ± 60 ± syst 1.7 

X PT to 0(p 4 ), "Fit II" in Ref. [16] 938 (fixed) 

XQSM (only L > 2 fm) 927 ± 4 ± syst 0.5 

XQSM (all data points) 931 ± 4 ± syst 2.3 



finite range regulator (FRR) method of Ref. [24] . There 
5 free parameters appear, which can well be constrained 
by lattice data thanks to the larger range of applicability 
of the approach, namely up to to 2 = 1 GcV 2 . In Ref. [24] 
different shapes of regulators were exploited and shown 
to yield practically the same results. Fig. 5b shows the 
5-parameter fit (using the dipole-type regulator) to the 
CP-PACS lattice data sets with (3 = 2.1 and 2.2 from 
Ref. [4]. In the xQSM-approach with one free parameter 
only, we were able to fit both data sets separately, see 
Figs. 2c and d. For sake of comparison we include the 
combined xQSM-fit in Fig. 5b and summarize the results 
in Table IV. 

As demonstrated in Fig. 5b and Table IV the results 
for Mjv(m^) from the FRR approach and the xQSM de- 
scribe the CP-PACS data [4] equally well. Having a closer 
look on the region to 2 < 0.3 GcV 2 we see that M^TOtt) 
from the xQSM starts to deviate more and more strongly 
from the fit of the FRR-approach with decreasing to 2 , 
though the curves remain very similar. The reason for 
the discrepancy (which apparently was of no relevance 
for to, 2 > 0.3 GeV 2 , see Fig. 5b) should be attributed 
to the fact, that Ma — Mm is kept finite in the FRR 
approach, but neglected in the xQSM. 

Thus one is lead to the conclusion, that the FRR 
method and the xQSM are completely consistent modulo 
1/N C corrections. 

At the physical point the difference between the values 
of Mm extrapolated by means of the FRR method and 
the xQSM is about 40 MeV, which again may give us a 
rough idea on the theoretical uncertainty due to neglect- 
ing the A-nucleon mass difference. Note that for this 
rough estimate we use the central value of the extrapo- 
lated Mm from the FRR approach, see Table IV, which 
has a substantially larger statistical uncertainty arising 
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TABLE IV: The value of the physical nucleoli mass Mjv as extrap- 
olated from CP-PACS lattice data [4] obtained from lattices of the 
sizes L = (2.2 — 3.1) fm with lattice spacings a = (0.09 — 0.13) fm 
covering the range 520 Mev < < 970 MeV. As guidelines for 
the extrapolation were used the finite range regulator approach [24] 



and the \QSM. Sec also the remarks 


in the caption to Tabic 


■ III. 


Method 


Mjv in McV 


xl ot 


finite range ("dipole") regulator [24] 


959 ± 116 ± syst 


0.4 


XQSM 


920 ± 3 ± syst 


1.2 



from fitting 5 parameters to the lattice data. In this re- 
spect the xQSM-fit is far more precise. 

To summarize, the comparison of the xQSM-based fits 
and those obtained using the first-principle approach in 
XPT [7, 16] and the FRR approach [24] leads us to the 
following conclusion. The systematic uncertainty of Mm 
due to neglecting the finite A-nucleon mass-splitting in 
the xQSM is effectively of the order of magnitude of 
50 MeV with the tendency to underestimate the nucleon 
mass. Noteworthy, a similar result - concerning both sign 
and order of magnitude - follows from a crude estimate 
within the model itself, see App. C. 

Recall that we did not consider isospin breaking ef- 
fects or electromagnetic corrections ~ O(10 MeV), see 
footnote 1. Including this we have to assign a systematic 
error to the xQSM-fits of about (<5M i v) S yst ~ 60 MeV. 
Taking into account a systematic error of this magnitude 
we observe that the extrapolated values in Tables II, III 
and IV are all consistent with the physical mass of the 
nucleon. 

We stress that we do not see any possibility to quan- 
tify the uncertainty of the xQSM-extrapolation of lattice 
data more quantitatively than that. 



VI. CONCLUSIONS 

The implicit dependence of the nucleon mass Mpj on 
the pion mass m w was studied in the large- N c limit in 
the framework of the chiral quark soliton model. The 
M^f(m v ) in the model exhibits a chiral behaviour and 
includes leading non-analytic terms which are consistent 
with the large- A^ formulation of QCD [35]. As was shown 
here, the model describes correctly also the heavy quark 
limit. The most remarkable observation we make here is 
that the model results for Mjv(m T ) well describe lattice 
data from full simulations [3-9] over the wide range of 
pion masses 0.1 GeV 2 < to 2 < 2.5 GeV 2 , provided one 
takes into account the generic overestimate of the nu- 
cleon mass in the soliton approach [38]. This is done by 
introducing an TOTr-indcpcndcnt subtraction constant, i.e. 
one single parameter to be fitted to the respective lattice 



data set. 

The good description of the lattice data on M^im^) 
can partly be understood as follows. In the xQSM, in the 
leading order of the large- N c limit, the A-nucleon mass- 
splitting A = M A - M N ~ 0{N~ 1 ) is neglected. That 
is a reasonable approximation when comparing to the 
present day lattice simulations where A 2 <C to 2 holds. 
However, the remarkable precision, to which the model 
describes the lattice data, remains a puzzle — to be clar- 
ified by further model studies. 

We observe that the values for the nucleon mass "ex- 
trapolated" from the lattice data [3-9] on the basis of 
results from the xQSM are in good agreement with ex- 
trapolations based on the first principle approach in xPT 
[7, 16] or the effective FRR approach of Ref. [24], and 
agree with the physical nucleon mass to within 5%. (But 
one has to keep in mind that at the same time — for the 
adopted handling of model parameters in the chiral limit 
- the lattice values of the pion decay constant at large 
TOtt <~ 1 GeV are underestimated by up to 40 %.) 

It is difficult to exactly quantify the theoretical uncer- 
tainty of this extrapolation due to the model-dependence. 
The main limitation for using the xQSM as a guideline 
for the chiral extrapolation of lattice data is due to the 
large- N c limit. There is no strict control within the model 
of the theoretical uncertainty introduced by neglecting 
the finite A-nucleon mass-splitting at the physical point, 
and we cannot quantify this and other uncertainties due 
to model-dependence quantitatively. This limits the use 
of the model as a quantitative effective tool for the ex- 
trapolation of lattice data. 

Still, the model may provide interesting qualitative in- 
sights — in particular in those cases when the matching 
between lattice results and xPT is difficult. A prominent 
example for that are (moments of) structure functions. 
In order to use the xQSM as a qualitative, but within 
its model accuracy reliable device, which is helpful for a 
comparison of lattice results to experimental data, fur- 
ther model studies are necessary — concerning the issue 
of handling model parameters in the chiral limit as well 
as addressing other observables [44] . 

From the model point of view the observations made 
in this work also contribute to a better understanding 
of the physics which underlies the chiral quark-soliton 
model, and may — that is our hope — stimulate further 
studies in this direction in this, and perhaps also other 
models. 
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APPENDIX A: PROPER-TIME 
REGULARIZATION 

The proper-time regularized versions of the integrals 
Ji and I2 appearing in Eqs. (2, 4) are given by 



h(m) 



du exp(— uM' ) 
t? (4^)2 



(Al) 



,'dit exp(— uM' ) 
L " ini] = I to (4tt)2 /d/3exp(n/3(l-/3K) 



where the m-dependence is hidden in M' = M + to. For 
a given and / w = 93 MeV fixed and due to Eqs. (4, 5) 
both the current quark mass and the cutoff are (implicit) 
functions of to t , i.e. to = m(m w ) and A cu t = Acut^ir); 
some selected values are shown in Table I. 

The expression for the soliton energy (9) in the proper- 
time regularization is given by 



£sol = N, 
R(u) 



E lcv +J2(R(E n ) - R(E no )) 

n 

00 

d u 



40F J u 3/2 
a; 



cxp(— UU) ) . 



(A2) 



APPENDIX B: THE PION-NUCLEON 
SIGMA-TERM a vN 

The pion-nuclcon sigma-term is an important quantity 
to learn about chiral symmetry breaking effects in the 
nucleon. The Feynman-Hellmann theorem [47] relates 
oVat to the slope of Mat (to) (with m = m q = m u = rrid 
neglecting isospin breaking effects) as follows 

dM N (m) 2 dM N (m v ) 
a wN =m — = r< — ^ (Bl) 



dm 



dmi 



where the second equality holds, strictly speaking, only 
for small m v . Eq. (Bl) offers a convenient way to learn 
about ovjv from lattice calculations of Mjv(to^), see e.g. 
[22]. Direct lattice calculations are more difficult but 
possible, see e.g. [48]. 

Apart from (Bl) one also can evaluate in the model 
directly the actual definition of <j„n as double commuta- 
tor of the strong interaction Hamiltonian with two axial 
isovector charges [49, 50], or exploit a sum rule for the 
twist-3 distribution function e(x) [51]. All three methods 
yield the same result in the xQSM [35]. 



By means of Eq. (Bl) we obtain in the present study 
in the proper-time regularization: g^n — 40 MeV. This 
agrees well with other xQSM calculations [30, 49, 50] 
performed in this regularization scheme. 

In the xQSM the pion-nucleon sigma-term is quadrat- 
ically UV-divergcnt, and as such particularly sensitive to 
details of regularization. In Ref. [35] a way was found 
to compute cr^N in a regularization scheme independent 
way with the result ovjv ~ 68 MeV. The prize to pay 
for the regularization scheme independence was the use 
of an approximation expected to work within an accu- 
racy of 0(30%) and well-justified by notions from the 
instanton vacuum model. 

On the basis of our results for Mjv(to^) we are now 
in a position to check the accuracy of this approxima- 
tion in practice. For that we note that the Feynman- 
Hellmann relation (Bl) allows to determine Mjv(to^) 
from (T,jj(m,) up to an integration constant. 

From the approximate (but regularization scheme in- 
dependent) result for oVA^m^) from Ref. [35] one obtains 
for Mjv(m 7r ) the result shown as dashed line in Fig. 1, 
where for convenience the integration constant is cho- 
sen such that both curves coincide at a central value of 
to^ = 1 GeV in the plot. We observe a good agreement 
of our exact result and the approximation of Ref. [35], 
see Fig. 1. 

Also our result for o^n agrees with the regulariza- 
tion scheme-independent result for u^n from Ref. [35] 
to within the expected accuracy of O(30%). 

Thus, our results confirm that the instanton model 
motivated approximation advocated in [35] works well. 
Note, however, that here we went far beyond the study 
of Ref. [35]. In this work we practically demonstrated 
that stable soliton solutions do exist also for large values 
of the pion mass, while in Ref. [35] this was presumed. 

The regularization scheme independent result of 
Ref. [35] is in good agreement with recent extractions 
indicating ct^n = (60 — 70) MeV [52] , which is substan- 
tially more sizeable than the value obtained from earlier 
analyses [53]. 

In this context it is worthwhile mentioning that the 
spectrum of exotic (" pentaquark" ) baryons which were 
predicted in the framework of the xQSM [54] and for 
which recently possible observations were reported - see 
[55] for recent overviews - could provide an independent 
mean to access information on the pion-nucleon sigma- 
term [56]. The presently available data on the exotic 
baryons favour a large value for ovat as found in [52]. 



APPENDIX C: SYSTEMATIC UNCERTAINTY 
OF THE ROTATING SOLITON APPROACH 

The neglect of A = Ma — Mjv in the leading order 
of the large- N c limit introduces a systematic uncertainty 
which is difficult to quantify. Here we roughly estimate 
this uncertainty on the basis of the soliton approach. 

In this approach the nucleon and the A-resonance are 
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just different rotational excitations of the same classic 
object, the soliton. A non-zero mass difference between 
the nucleon and the A-resonance arises due to consider- 
ing a particular class of ("rotational") l/iV c -corrections 
[25, 26]. The mass of an SU(2)-baryon M B (B = N or 
A in the real world with N c = 3 colours) with spin Sb is 
given by 

M B = E sol + M 2 + Sb(S b + + (CI) 

with E so \ as defined in Eq. (9), M 2 denoting as in Eq. (12) 
the 0(N®) correction to the baryon mass which is the 
same for the nucleon and the A-resonance, and the dots 



representing higher order l/7V c -corrections. 

Then we obtain as a first correction to Mm in Eq. (22) 



M JV (m w ,A)=M JV (m*,0) + -Y^ + ... (C2) 

Focusing on the linear order correction in A and neglect- 
ing higher orders, we see that by neglecting A one un- 
derestimates the nucleon mass by about 70 MeV at the 
physical value of the pion mass. This is in good agree- 
ment with the systematic uncertainty roughly estimated 
in Sec. IV. 
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